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In this paper we develop a white noise framework for the study 
of stochastic partial differential equations driven by a d-parameter 
(pure jump) Levy white noise. As an example we use this theory to 
solve the stochastic Poisson equation with respect to Levy white noise 
for any dimension d. The solution is a stochastic distribution process 
given explicitly. We also show that if d < 3, then this solution can 
be represented as a classical random field in L^(/x), where n is the 
probability law of the Levy process. The starting point of our theory 
is a chaos expansion in terms of generalized Charlier polynomials. 
Based on this expansion we define Kondratiev spaces and the Levy 
Hermite transform. 



1. Introduction. White noise analysis has become a subject of much 
current interest. This theory was first treated by Hida [14] and extensively 
studied in many other works. See [16] and the references therein. These 
investigations are based on the concept of a Gaussian measure and the asso- 
ciated expansion into Hermite polynomials. Later on an extension of white 
noise theory to non-Gaussian analysis was established in [4] and developed 
further by Kondratiev, Da Silva, Streit and Us [24] and Kondratiev, Da Silva 
and Streit [23]. The main tool of this theory is a biorthogonal decomposi- 
tion, which extends the Wiener-Ito chaos expansion. White noise analysis 
has been used in a broad range of applications. This approach was originally 
applied in quantum physics. See, for example, [3] or [2]. Subsequently, new 
applications have been found in stochastic (partial) differential equations 
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[17]. See also [21] and [6] to mention a few. More recently, the theory has 
been applied to finance [1]. See [18] and [12] for the fractional Brownian 
motion case and [10] and [34] in the non-Gaussian case. 

The object of this paper is to provide a white noise framework, based 
on results in [28, 10, 34] and [17], to study SPDEs driven by (pure jump) 
Levy processes. We apply this theory to solve the stochastic Poisson equation 
driven by a d-parameter ( pure jump) Levy white noise. That is, consider the 
following model for the temperature U{x) at point x in a bounded domain 
D in M"^. Suppose that the temperature at the boundary dD of D is kept 
equal to zero and that there is a random heat source in D modeled by Levy 
white noise V{x) =V{xi, . . . , x^)- Then U is described by the equation 

, . AU{x) = -Vix), x = {xi,...,Xd)£D 

^ ■ ^ U{x) = 0, xedD. 

It is natural to guess that the solution must be 

(1.2) U{x) = U{x,uj)= f G{x,y)dr]{y), 

Jd 

where G{x,y) is the classical Green function for D and the integral on 
the right-hand side is a multiparameter Ito integral with respect to the 
d-parameter Levy process rj{x). But the integral on the right-hand side of 

(1.2) only makes sense if G{x, •) is square integrable in D with respect to 
the Lebesgue measure. The latter is true if and only if the dimension d is 
chosen lower than 4. Despite this difficulty we will show the existence of a 
unique explicit solution 

X U{x, •) G (5)-i, 

where (5)_i is a suitable space of stochastic distributions, called the Kon- 
dratiev space. 

The stochastic Poisson equation (1.1) was discussed by Walsh [39] in the 
case of Brownian white noise W. He proved that there exists for all d a 
Sobolev space H~''^{D) and an i7~"(D)-valued stochastic process 

U = U{uj):^^H-''{D) 

such that (1.1) holds in the sense of distributions, for example, 

([/(•, cj), A(/.) = -(W^(-,a;),0) a.s. for all (/> G i7~"(L>). 

The solution of Walsh is given explicitly by 

(1.3) {U,<P)= [ I G{x,y)<i^{x)dxdB{y)- <j>&H-^{D). 
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The system (1.1) was also studied in [17] in the Gaussian case. There the 
solution U{x), which takes values in the Kondratiev space, can be described 
by its action on test functions / € (5)i, 

(1.4) {U{x)J)=[ G{x,y){W{y),f)dy; /g(5)i. 

If we compare (1.3) and (1.4) we find that the Walsh solution takes x- 
averages for almost all uj, whereas the last one takes w-averages for all x. 

Our solution is an extension of (1.4) to Levy processes. The approach we 
use to solve (1.1) is based on a chaos expansion in terms of generalized Char- 
lier polynomials (cf. [28]) and on concepts developed in [17, 10] and [34]. Our 
method, which can be applied to other classes of SPDEs, has the advantage 
that SPDEs can be interpreted in the usual strong sense with respect to 
time and space. There is no need for a weak distribution interpretation with 
respect to time and space. Furthermore, the Walsh construction reveals the 
disadvantage of defining a multiplication of (Sobolev or Schwartz) distribu- 
tions, if one considers SPDEs, where the noise is involved multiplicatively. 
However, on the Kondratiev space (5)_i we can define a multiplication, the 
Levy Wick product. This gives a natural interpretation of SPDEs, where 
the noise or other terms appear multiplicatively. Furthermore, in some cases 
solutions can be explicitly obtained in terms of the Wick product. See [17]. 

The general machinery, developed in this paper, is of independent interest 
and we are convinced that it serves a useful tool for the study of a large class 
of stochastic partial differential equations driven by Levy space-time white 
noise. 

Finally, let us mention that there has recently been an increasing interest 
in solving SPDEs driven by d-parameter Levy processes. We refer to [5, 29] 
and the references therein. 

We shall give an overview of the paper. In Section 2 we introduce a white 
noise framework for the study of SPDEs driven by d-parameter Levy pro- 
cesses. The starting point of our theory is a chaos expansion in terms of 
generalized Charlier polynomials. Based on this expansion we define Kon- 
dratiev spaces, the Wick product and the d-parameter Levy white noise. 
Further, we give the definition of the Levy Hermite transform and state a 
characterization theorem for the Kondratiev space (5)_i. In Section 3 we 
use the tools developed in Section 2 to apply it to solve the stochastic Pois- 
son equation driven by a d-parameter Levy white noise. Finally, in Section 4 
we discuss the solution and its properties. In particular, we show that if 
d <3, the solution can be represented as a classical random field in L^(/x), 
where /i is the underlying probability law of the Levy process. 

2. Framework. In this section we give the general framework to be used 
later. The starting point for our discussion are white noise concepts for Levy 
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processes, developed in [10, 34] and [28]. Actually, we empasize the use of 
multidimensional structures, that is, the white noise we intend to consider is 
indexed by a multidimensional parameter set. Our presentation and notation 
will follow that of [17] closely, where Gaussian white noise theory is treated. 
For more information about white noise theory we refer to [16, 26] and [32]. 

2.1. A white noise construction of Poisson random measures associated 
with a Levy process. In this paper we confine ourselves to (d-parameter) 
pure jump Levy processes without drift. 

A pure jump Levy process r]{t) on M with no drift is a process with in- 
dependent and stationary increments, continuous in probability and with 
no Brownian motion part. The characteristic function of such a process is 
given by the Levy-Khintchine formula in terms of the Levy measure u of 
the Levy process, that is, in terms of a measure v on Mq ■.= M — {0}, that 
integrates the function 1 A z^. Hence, driftless pure jump Levy processes 
can be characterized as Levy processes with characteristic triplet (0,0,1/). 
For general information about Levy processes see [8] and [36]. In general, 
such processes do not possess the chaotic representation property, but they 
admit a chaos representation with respect to Poisson random measures (see 
[19]). Therefore, we aim at viewing these processes as elements of a certain 
Poisson space. In this framework we will give a white noise construction of 
Poisson random measures and, since our emphasis lies on processes indexed 
by multidimensional sets, we will define d-parameter (pure jump) Levy pro- 
cesses. Further, we prove a chaos expansion in terms of generalized Charlier 
polynomials. 

A usual starting point in white noise analysis is the application of the 
Bochner-Minlos theorem to prove the existence of a probability measure 
on the space of tempered distributions 5'(]R'^). However, it turns out that 
5'(IR'^) is not the most appropriate for dealing with Levy processes since 
this choice would require restrictive conditions to be imposed on the Levy 
measure. This circumstance comes from the fact that the Levy measure has 
a singularity at zero. Therefore, we use the construction of a nuclear algebra 
S{X), which is more tractable for our purpose. In fact, the space S{X) is 
a variant of the Schwartz space on X = M'^ x Mq, more precisely S{X) is a 
subspace of the Schwartz space modulo a certain sub^pace depending on the 
Levy measure. Let us first give the construction of S{X) (cf. [28]). 

In the following let {^n}n>o be the complete orthogonal system of L^(M), 
consisting of the Hermite functions. Then the (countably Hilbertian) nuclear 
topology of the Schwartz space 5(M'') is induced by the compatible system 
of norms 

(2.1) yf^:= ^(l + a)2^(v,,Ui2(K.), 7G<, 
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where := Ilti U and (1 + a)^^ := nti(l + "i)^^' for a = (ai, . . . , a^) G 
N*^ and 7 = (71, . . . , 7^;) S Nq. Now let us take a numbering of the norms in 
(2.1) 

) say II • ll-y. , and define the norms \\ip\\p — X]f=i llfW^i > P G N. Then || • ||p, 
p € N are increasing pre-Hilbertian norms on It is weU known that 

the norms || • ||p are equivalent to the norms || • ||g,oo) given by 

(2.2) ||9?||5,oo:= sup sup |(l + |z|'=)a^(/j(z)|, g S Nq, 

0<k,\'r\<q zeR'i 

where d'^ip = q^jtz^'P for 7 = (7i, • • • , Id) G ^0 with I7I := 71 + • • • + 7^. 

1 d 

We mention the following important property of the norms || • ||p (see [20]): 
For all p G N there exists a constant Mp > such that for all (p,ip G S{M.'^), 

(2.3) ||(/?V||p < ^'^plbllpllV'llp- 
We then define the space S{X) by 

(2.4) S{X) := € 5(M'^+i) : ip{zi,. ..,Zd,0) = (^^'/') ■ • ■ > ^d, 0) = o| 

It follows that S{X) is a closed subspace of 5(]R'^"''^). Thus, S{X) is a (count- 
ably Hilbertian) nuclear space with respect to the restriction of the norms 
II • lip. Moreover, it is a nuclear algebra, that is, S{X) is, in addition, a topo- 
logical algebra with respect to the multiplication of functions. In the sequel 
we denote by A^'^ the Lebesgue measure on R*^ and by a Levy measure of 
on Mq- We set vr = A^"^ x u. We need the following result. 



Lemma 2.1. There exists an element denoted by 1(8) in S'{X) such 
that 

(1® / 0(y)7r(dy) 

Jx 

for all (j) G S{X), where (1(8) = (1(8 z^)(0) is the action of 1(8 on (j). 
The notation v shall indicate that u is the Radon-Nikodym derivative of v 
in a generalized sense. 

Proof. Without loss of generality we consider the case d=\. Set L((^) = 
ip{z)Tr{dz). Let ipn, (p € S{X) with (pn^ (p in S{X). By Taylor's formula 
we have for ip £ S{X) that 

93(2;, z) = ip{xO) + (^^V') {^^ 0)2; + ^ (^■^'P^ (x, i)z'^ 
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for a point ^ between and z. We assume w.l.o.g. that the measure v vanishes 
outside of [-1,0) U (0, 1]. Therefore, it follows by (2.2) that 

\L{^pn - ^)\< [ [ \(pn{x,z) - (f{x,z)\i/{dz)X{dx) 
JrJ-1 



< 



1 (1 + |xp + \z\^)\ipn{x,z) - ip{x,z) 
-1 



z^ 



(1 + \x 



-v{dz)\[dx) 



< \\v>n - ^h.oo [ /-I /i lo^ Hdx) [ z^v{dz) 
■ Jr(1+|x|^J J -I 

— > for n — > CO. 
So the linear functional L is continuous on S{X). □ 

Next define the space 

(2.5) K:={<P(^S{X):U\\l2^^)=Q}. 

By the same arguments as in the proof of Lemma 2.1 it can be shown that 
A^TT is a closed subspace of S{X). Furthermore, one checks that it is a closed 
ideal of S{X). Now we introduce the space 5(X), which we use to construct 
the white noise measure. 

Definition 2.2. We define the space S{X) as follows, 

(2.6) S{X)=S{X)/M^. 

The space S{X) is a (countably Hilbertian) nuclear algebra with the com- 
patible system of norms 

(2.7) ||^||p,^:= inf ll^ + V'llp, j> G N. 

See page 72 in [13]. Further, let S\X) denote the topological dual of S{X). 
We obtain the following corollary to Lemma 2.1: 

Corollary 2.3. The functional L{(f)) := J-^ (j)(z)7r(dz) satisfies the in- 
equality 

\L{$)\ < MpU\k^ 
for all p > po, which yields the continuity of the functional L on S{X). 
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Theorem 2.4. There exists a unique probability measure fi on the Borel 
sets of S'{X) with the following Poissonian characteristic functional with 
intensity vr such that for all (j) G S{X): 

(2.8) / e^<'^''^> dfi{u;) = exp f / (e^"^ - 1) dn) , 

JS'(X) \Jx J 

where {uj, (j)) = uj{(j)) is the action of € S\X) on (p ^ S{X). Moreover, there 
exists a Po £ N such that 1 ® z> G 5_pq(X) and a natural number qq > pQ 
such that the imbedding operator Sqg{X) ^ Spf^{X) is Hilbert-Schmidt and 
^(5_g,j(X)) = 1. The space Sp{X) denotes the completion of S{X) with re- 
spect to II • llp^TT O'nd S-p{X) is the corresponding dual with norm \\ ■ ||-p,7r. 

Proof. Since |e*^ ~ 1| 1^ \z\, the result follows from Corollary 2.3 and 
Bochners theorem for conuclear spaces [13]. □ 

We call the probability measure fi on 0, = S'{X) in Theorem 2.4 Levy 
white noise probability measure. It turns out that this measure satisfies the 
first condition of analyticity in the following sense (see [23] ) . 

Lemma 2.5. The Levy white noise measure satisfies the first condition 
of analyticity, that is, there exists e > and a po such that 



I exp(e||a;||_pQ) dniui) < oo. 
JS'(X) 



IS'iX) 

Proof. The proof follows the argument of Lemma 3 in [38]. Introduce 
the moment functions of /x, which by a criterion of Cramer [9] can be ex- 
pressed by 



M„(0) := / {u;,cPrd^,{co) = ^L{t(t>) 



t=o 



for every <j) E S{X), n G N. Define the set 



A^:= |(ai,...,afc) eN^^ 
I i=i 



ai = n , 



Then we obtain the following expression for 

(2.) M„w^x:^Efi^«. 

fc=i aeA,^; i=i ^ 
We get for the number pQ in Theorem 2.4 that 

1(1 (g) z>,0)| < ||1 (g)z>||_po,^||(/)||po,^ < oo. 
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Next relation (2.3) implies that for all p G N there exists a constant Mp > 
such that for all (pjip £ '5(X), 

(2.10) \\mp,n<Mp\mpMp,^- 

Thus, we get that 

|(l®z>,(/."^)| < \\l^O\\_p,,^{Mp,pMZ^, 
if we choose Mp^ > 1. So we deduce from (2.9) that 

fc=l asA^ i=l ^ 

= F„(||l0z>||„p„,^)C,"J|</.||^„,,, 

where Fn{x) is the nth moment of the Poisson distribution with intensity 
X and where Cp^ is a constant. Further, it is known that for a Poisson 
distribution with intensity x= ||1 t'||-pQ.7r, there exists a constant Cx such 
that for ah n G N, 

|F„(||l0z>||„p„,OI<n!C^i^,||_^^^^. 
Therefore, we get for a C > that 

|M„(0)|<n!C"M;,,. 
The claimed result follows from Lemma 3 in [23]. □ 

Further, consider the function a defined by a{cj)) = log(l + ip) modA/'Tr for 
4> = ip with (p{x) > — 1. Note that a is holomorphic at zero and invertible. 
With the help of Lemma 2.5, it can be shown just as in [28] that there exist 
symmetric kernels C„(w) such that for all (j) in an open neighborhood of zero 
in S{X), 

where c/*'^" G S{X)®'^ . The symbol S{X)'^'^ denotes the nth completed sym- 
metric tensor product of S(X) with itself. The elements of this space can 
be seen as functions / G 5(X"') modulo 7V^x»i such that / = /(xi, . . . ,x„) 
is symmetric with respect to the variables xi, . . . ,x„ G X. From (2.11) we 
conclude that the Cn are generalized Charlier polynomials (see [23]). We 
have that 

(2.12) {(^(a;), c^^^^) : ^^^^ G 5(X)^", n G No} 
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is a total set in -L^(^). Furthermore, for all n, m, (p^"'^ e and i/j^^^ G 

the orthogonality relation 

(2.13) / (C„M,(^("))(C^(a;),V'("^))d^(a;) = 5„,^n!(^W,V'W)^2(^„) 
holds. See [28]. 

Remark 2.6. It can be easily seen from (2.13) and the construction of 
S{X) that the Levy white noise measure fi is nondegenerate in the following 
sense (see [23]): Let F be a continuous polynomial, that is, F is of the 
form F{lo) = E^=i(^®", f^^"^) for u € S'{X), iV G Nq with e 5c(X)®" 
[complexification of 5(X)®"]. If F = n-a.e., then -F(u;) = for all uj G 
S'{X). We mention that this property is essential for the construction of 
certain test function and distribution spaces (see [23, 28]). 

Next, for functions / : — > M define the symmetrization (f)^ of / by 

(2.14) (/)^(xi, . ..,Xr,):=-J2 /(^-i , • • • , ^.n) 

n 

for all permutations a of {1, . . . , n}. Then a function / : X^ ^ R is symmet- 
ric, if and only if f = f. Denote by L^(X"', vr^") the space of all symmet- 
ric functions on X^, which are square integrable with respect to vr^". Let 
fn G L2(X",7r^"). Since S{X) is dense in L'^{X,tt) (cf. [28]), we can choose a 
sequence fi'^ in 5(X)®" with fH^ fn in L2(X",7r^"). Then (2.13) implies 
the existence of a well defined {Cn{oj)-, fn) such that 

(2.15) (a(a7),/„) = lim(a(u;),/«) in vr><"). 
Since Ci(a;) = cj — 1 (g) t' for all G 5(X) (see [28]), we get 

(2.16) / (u;-l®z>,/)2(iM^) = ||/||i,(^). 

Further, if we define for Borelian Ai C M*^, A2 C Mq with 7r(Ai x A2) < c3o 
the random measures 

iV(Ai, A2) := (a;,XAixA2) and 

(2.17) 

iV(Ai,A2) := (i^ - l®z>,XA,xA2), 

we see from their characteristic functions that is a Poisson random mea- 
sure and is the corresponding compensated Poisson random measure. The 
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compensator of iV(Ai,A2) is given by vr. Therefore, it is natural to define 
the stochastic integral of cp £ L'^i'^) with respect to N by 

(2.18) j (l){x, z)N{dx, dz) := {oj 
In particular, if we define 

■■= / X[o,xi]x-x[o,xa]ix) ■ zN{dx,dz) 

(2.19) ^ 

for X = (xi , . . . , Xd) G W^, 

where [0,Xj] is interpreted as [xj,0], if Xj < and where the Levy measure 
u is assumed to integrate z'^, then fj{x) has a version r/(x), which is cadlag 
in each component Xj. This follows with the help of (2.13). We call r/(x) 
d-parameter Levy process or space-time Levy process. 

We conclude this section with a chaos expansion result in terms of the gen- 
eralized Charlier polynomials C„. The result is a consequence of (2.12) and (2.13). 

Theorem 2.7. If F £ L'^^j,), then there exists a unique sequence fn G 
L2(X") such that 

(2.20) F{u;) = J2{Cn{u;)Jn). 

n>0 

Moreover, we have the isometry 

(2.21) \\F\\W) = Y.^'-\\fn\\hix^y 

n>0 

2.2. Chaos expansion, Kondratiev spaces {S)p, {S)-p and Levy white noise. 
First we reformulate the chaos expansion of Theorem 2.7. Then we use the 
new expansion to define a Wick product on spaces of stochastic test functions 
and stochastic distributions. The definitions and results here are analogous 
to the one-parameter case, which is treated in [10] and [34]. 

From now on we suppose that our Levy measure u satisfies the condition 
of [30] , namely, that for every e > there exists a A > such that 

(2.22) / exp{X\z\)u{dz) <oo. 

JR\(-e,e) 

This implies that our Levy measure has finite moments of all orders > 2. 

For later use we introduce multi-indices of arbitrary length. To simplify 
the notation, we regard multi-indices as elements of the space (Nq )c of all 
sequences a = {ai,a2, ■ ■ ■) with elements Ui G Nq and with compact support, 
that is, with only finitely many 7^ 0. We define 

J = {K)c. 
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Further, we set Index(a) = max{i : 7^ 0} and |a| = J2i for J . 

Next we consider two famihes of orthogonal polynomials. We use these 
polynomials to reformulate the chaos expansion of Theorem 2.7. First let 
{Cfc}A:>i be the Hermite functions, just as in Section 2.1. Now choose a 
bijective map 

/i : N"* ^ N. 

Define the function Cfc (xi , . . . , x^) = (xd) , if /c = , . . . , z^) for G 

N. Then {Cfc}fc>i constitutes an orthonormal basis of L^(R'^). 

Further, let {lm}m>o be the orthogonalization of {1, z, z'^ , . . .} with re- 
spect to the innerproduct of L'^{q), where g{dz) = z'^v{dz). Then define the 
polynomials 

(2.23) p^[z) = — \ z-lm-i{z). 

IKm-l||L2(p) 

The polynomials pm form a complete orthonormal system in L'^{v) (see 
[34]). We shall mention that we could also use any orthonormal basis in 
S{X) C L'^{v) for d = instead of the polynomials pm- In this case the 
integrability condition (2.22) reduces to the requirement of the existence of 
the second moment with respect to v. The choice of the polynomials pm 
serves to ease notation. 

Next define the bijective map 

(2.24) z:NxN^N; (i, j) ^ j + (i + j - 2)(i + j - l)/2. 

Note that z{i,j) gives the "Cantor diagonalization" of N x N. 
Then, if fc = z(i,j) for i,j G N, let 

6k{x,z) = Ci{x)pj{z). 

Further, assume Index(Q) = j and \a\ = m for J and identify the func- 
tion 5*^" as 

(5^"((xi,zi),...,(xm,Zm)) 

= Ci®---05f"^((xi,zi),...,(x^,z^)) 

(2.25) 

= (5i(xi,2;i)---5i(x«i,z„J 

• ■ • (^j (Xq-^-i hOj-i+l 5 •^aiH hoj-i+l) ' ' ' {•^m i Zrn) ; 

where the terms with zero-components Qj are set equal to 1 in the product 
{6f = l). 

Finally, we define the symmetrized tensor product of the 5k s, denoted by 
as 

6^''{{xi,Zi),...,{Xm,Zm)) 

(2.26) =(5^-)^((xi,zi),...,(x^,z™)) 

= ^r^® • ■ • 0<5f ((Xi, Zi), . . . , (Xm, Zm)). 
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For J define 

(2.27) K„H:=(q,|(a;),0, 
where we let Kq{uj) = 1. For example, if a = with 

(2.28) .'(i) = {J; '^^''i>i. 

we obtain 

(2.29) K,^{u:) = (u;,5®^') = {lo,5i) = {u;,Qix)p,{z)), 
ifl = z{i,j). 

By Theorem 2.7 any sequence of functions fm G L'^ {tt^"^) , m = 0,1,2, ... , 
such that X]m>i "^'Il/m|li2(7rxm') < oo defines a random variable F G L'^in) 
by -F(a;) = J2m>o{Cm{^^) , fm) ■ Since each is contained in the closure of 

the linear span of the orthogonal family {6'^°'}\a\=m in L'^{'7t^''^), we get for 
all m > 1 the representation 

(2.30) /m= E C"-^®" 

\a\=m 

in L^(7r^™') for G M. Hence, we can restate Theorem 2.7 as follows. 

Theorem 2.8. The family {Ka}aej constitutes an orthogonal basis for 
L'^ifJ') with norm expression 

(2.31) Il-^'^"alli2(;,) = Oil := ai\a2\ 

for a = (ai, a2, . . . ) € J". Thus, every F € L^ifJ^) has the unique representa- 
tion 

(2.32) F=Y,CaKa, 

where Co G K for all a and where we set cq = E[F]. 
Moreover, we have the isometry 

(2.33) = alcl 

Example 2.9. (i) Choose F{uj) =r/(x), the (i-parameter Levy process. 
Thenr?(x) =/[o^^^]^...^[o^^^]xiR^ zN{dx,dz) = (w,X[o,xi]x...x[o,xd(^)-^> a.e. and 
it follows by (2.29) that 



m ■■■ Ck{xi,...,Xd 
k>l ^0 ^0 



)dxi--- dxd- K^,(k,i), 
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where m = \\x\\i2(^,yy 

(ii) Let Ai C M™, A2 CMq with 7r(Ai x < 00. Set fi{x,z) = XAixAiix, z). 
Then by (2.29) and (2.30) we get for F = iV(Ai, A2) = (w, /i) 



(2.35) N{t,A)= I I C 



'k{x)pTri{z)v{dz) dx ■ K^^(k,m). 



Next we define various generahzed function spaces that relate to L?{lj) in 
a natural way. These spaces turn out to be a useful tool to study stochastic 
partial differential equations. Our spaces are Levy versions of the Kondratiev 
spaces, which were originally introduced in [22]. See also [4] and [25] in the 
context of Gaussian analysis. The one-parameter case with respect to the 
Levy white noise measure /i can be found in [10] and [34]. The extension to 
multidimensional parameter sets is analogous. 

Definition 2.10. (i) The stochastic test function spaces. Let < p < 1. 
For an expansion / = CaKa € L'^ifJ-) define the norm 

(2.36) \\f\\l,:=J2{a\)'+^clm''' 

for k G No, where (2N)'=° = (2 • l)''"i (2 • 2)'="2 • • • (2 • m)'="- , if Index(a) = m. 
Let 

(5)m:={/:||/||m.<oo} 

and define 

(2.37) {S)p:= fl (5),,fc, 

fcgNo 

endowed with the projective topology. 

(ii) The stochastic distribution spaces. Let < p < 1. In the same manner, 
define for a formal expansion F = ^^ej ^aKa the norms 

(2.38) := E (a!)^-''c2 (2N)-'=", k € Nq. 



Set 



and define 



(5)_p,_fc:={F: |lF|[_p,„;^.<oo} 



(2.39) {S)-p:= U ('5)-p,-fc, 

fceNo 

equipped with the inductive topology. 
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We can regard (<S)_p as the dual of {S)p by the action 

(2.40) {FJ)=J2baCaal 

for F = Y^aejbaKa G ('5)--p and / = Y^aej^aKa G {S)p. Note that for 
general < p < 1 we have 

(2.41) (5)i c (5)p c (5)o c l2(p) c (5)_o c (5)_p c (5)_i. 

The space (5) := (5)o, respectively, (5)* := (5)_05 is a Levy version of the 
Hida test function space, respectively, Hida stochastic distribution space. 
For more information about these or related spaces in the Gaussian and 
Poissonian case we refer to [16] and [17]. 

One of the remarkable properties of the space (<S)* is that it accomodates 
the (d-parameter) Levy white noise. See [10]. 

Definition 2.11. The (d-parameter) Levy white noise V{x) of the Levy 
process r]{x) (with m = ||^||i2(,^)) is defined by the formal expansion 

(2.42) V{x)=mY,Ckix)K^.^,,r^, 

k>l 

where Ck{x) is defined by Hermite functions, z{i,j) is the map in (2.24) and 
where € c7 is defined as in (2.28). 

Remark 2.12. (i) Because of the uniform boundedness of the Hermite 
functions (see, e.g., [37]) the Levy white noise f?(x) takes values in (5)* for 
all X. Further, it follows from (2.34) that 

(2.43) ^\ lix) =V{x) in (5)*. 

This justifies the name white noise for ^(x). 

(ii) Just as in [34] the {d-parameter) white noise N{x,z) of the Poisson 
random measure N{dx, dz) can be defined by 

(2.44) N{x,z)= J2 Ck{x)pm{z)-K.(k,^), 

fc,m>l 

where Pm{z) are the polynomials from (2.23). We have that N{x,z) is con- 
tained in (5)*7r-a.e. The relation (2.35) admits the interpretation of N{x,z) 
as a Radon-Nikodym derivative, that is, (formally) 
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The last relation entitles us to call N{x, z) white noise. 
Moreover, ?7(x) is related to N{x,z) by 

(2.46) ??(x)= / zN{x,z)u{dz). 

Jr 

The relation above is given in terms of a Bochner integral with respect to i' 
(see [34]). 

2.3. Wick product and Hermite transform. In this section we define a 
{stochastic) Wick product on the space (5)_i with respect to the Levy white 
noise measure /i. Then we give the definition of the Hermite transform and 
apply it to establish a characterization theorem for the space (5)_i. 

The Wick product was first introduced by Wick [40] and used as a renor- 
malization technique in quantum field theory. Later on a (stochastic) Wick 
product was considered by Hida and Ikeda [15]. This subject both in math- 
ematical physics and probability theory is comprehensively treated in Do- 
broshin and Minlos [11]. Today the Wick product provides a useful concept 
for a variety of applications, for example, it is important in the study of 
stochastic ordinary or partial differential equations (see, e.g., [17]). 

The next definition is a d-parameter version of Definition 3.11 in [10]. 

Definition 2.13. The Levy Wick product FoG of two elements 
F=Y, aaKa, G=J2 bpKp e (5)_i with aa,bi3 gM 

is defined by 

(2.47) FoG= aabpK^+f3. 

Remark 2.14. Let /n = E|a|=nCa5®" e L2(7r><") and 5™, = E|/3|=m x 
^ £2(-^xm) according to (2.30). Then we have 

fn^gm = E E = E E ^-hS^' 

\a\=n\f^\=m \'y\=n+ma+f3=-/ 

in L2(7r^ ("+'")). Hence, 

(2.48) {Cn{uj),fn)o{Cm{(^),gm) = {Cn+m{uj) , fn§>gm) ■ 

Remark 2.15. A remarkable property of the Wick product is that it is 
implicitly contained in the Ito-Skorohod integrals. The reason for this fact 
is that if Y{t) =Y{t,uj) is Skorohod integrable, then (see [10]) 

(2.49) rY{t)6i]{t)= r Y{t)o'h{t)dt. 
Jo Jo 
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The left-hand side denotes the Skorohod integral of Y{t) and the integral on 
the right-hand side is the Bochner integral on (5)*. The Skorohod integral 
extends the Ito integral in the sense that both integrals coincide, if Y{t,uj) 
is adapted, that is, we then have 

(2.50) rY{t)5r,{t)= ^ Y{t)dr,{t). 
Jo Jo 

Note that a version of (2.49) holds for the white noise N{t, x), too (see [34]). 
The extension to the d-parameter case is given in [28]. 

Remark 2.16. It is important to note that the spaces (5)i, (5)_i and (5), 
(5)* form topological algebras with respect to the Levy Wick product o (for 
an analogous proof see [35] and [17]). For more information about the Wick 
product and Skorohod integration in the Poissonian and Gaussian case see, 
for example, [16, 17] and [31]. 

The Hermite transform, which appeared first in Lindstr0m, 0ksendal and 
Ub0e [27], gives the interpretation of (5)-i in terms of elements in the al- 
gebra of power series in infinitely many complex variables. This transform 
has been applied in many different directions in the Gaussian and Poisso- 
nian case (see, e.g., [17]). Its definition for (d-parameter) Levy processes is 
analogous. 

Definition 2.17. Let F = Y,aej(^aKa e {■S)~i with G M. Then the 
Levy Hermite transform of F, denoted by HF, is defined by 

(2.51) nF{z) = J2 G C, 

if convergent, where z = (zi, Z2, . . . ) G (the set of all sequences of complex 
numbers) and 



n ) 



(2.52) z" = z'^'zr---z'^ 
if a = {ai,a2, . . . ) G J", where Zj = 1. 

Example 2.18. We want to determine the Hermite transform of the 
d-parameter Levy white noise Vix). Since Vix) = mJ2k>iCkix)K^^(k,i) , we 
get 

(2.53) n{v){x, z)=mJ2 Ck{x) ■ 2.(^,1) , 

fc>i 

which is convergent for all z G (C^)c (the set of all finite sequences in C^). 
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One of the useful properties of the Hermite transform is that it converts 
the Wick product into ordinary (complex) products. 

Proposition 2.19. IfF, Gg (5)_i, then 

(2.54) niFoG){z)=niF)iz)-niG)iz) 
for all z such that 7i{F){z) and 7i{G){z) exist. 

Proof. The proof is an immediate consequence of Definition 2.13. □ 

In the following we define for < R,q < oo the infinite-dimensional neigh- 
borhoods Kg{R) in by 

(2.55) ir,(i?) = |(6,6, • • • ) G : E in'mr <r']- 

By the same proof as in the Gaussian case (see Theorem 2.6.11 in [17]) 
we deduce the following characterization theorem for the space (5)_i. 

Theorem 2.20. (i) If F = J2aej(^aKa G then there are q, Mq < 

oo such that 

(2.56) \nF{z)\<Y, |a,||z"|<MjE(2N)''"k"p) ^ 



for all z S 

In particular, TiF is a hounded analytic function on Kq{R) for all R < 
oo. 

(ii) Conversely, assume that g{z) = X^aej ^o-^" power series of z ^ 
such that there exist q < oo, 6 > with g{z) is absolutely convergent 
and hounded on Kg{6) then there exists a unique G € (5)_i such that TCG = 
g, namely, 

(2.57) G=Y.b^Ka. 

3. Application: the stochastic Poisson equation driven by space-time Levy 
white noise. Let us illustrate how the framework, developed in Section 2, 
can be applied to solve the stochastic Poisson equation 

, . AUix) = -V{x), xeD, 

^ ^ Uix) = 0, xedD, 

where A = X)fc=i is the Laplace operator in M'^, D is a bounded do- 

main with regular boundary (see, e.g.. Chapter 9 in [33]) and where Vix) = 
mJ2k>iCk{x)K^^(k,i) is the d-parameter Levy white noise (Definition 2.11). 
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As mentioned in the Introduction, the model (3.1) gives a description 
of the temperature U{x) in the region D under the assumption that the 
temperature at the boundary is kept equal to zero and that there is a white 
noise heat source in D. 

Note that AU{x) in (3.1) is defined in the sense of the topology on (5)_i. 

Now we aim at converting the system (3.1) into a deterministic partial dif- 
ferential equation with complex coefficients by applying the Hermite trans- 
form (2.51) to both sides of (3.1). Then we try to solve the resulting PDE, 
and we take the inverse Hermite transform of the solution, if existent, to 
obtain a solution of the original equation. Before we proceed to realize our 
strategy, we need the following result. 

Lemma 3.1. Suppose X and F are functions from D in (3.1) to (5)_i 
such that 

(3.2) AnX{x,z)=nF{x,z) 

for all (x, z) E Z) X Kq{6) for some q < oo, 5 > 0. 

Furthermore, assume for all j that ■^j7iF{x, z) is bounded on D x Kq{5), 

continuous with respect to x & D for each z G Kq{6) and analytic with respect 
to ze Kq{5) for all xeD. 
Then 

(3.3) AX{x) = F{x) for all xeD. 

Proof. Use repeatedly the same proof of Lemma 2.8.4 in [17] in the 
case of higher-order derivatives. □ 

Now, we take the Hermite transform of (3.1) and we get 

/3_4N Au{x,z) = -n{'n){x,z), xeD, 

u{x, z) = 0, X G dD, 

where u = TiU and 7i{'n){x,z) = rnJ2k>iCk{x) ■ Zz{k,i) for z G (C^)c [see 
(2.53)]. By comparing the real and imaginary parts of equation (3.4), one 
checks that 

(3.5) u{x,z)= [ G{x,y)-n{il){y,z)dy, 

where G{x,y) is the classical Green function of D with G = outside of D 
(see, e.g., Chapter 9 in [33]). Since G{x, •) G L-^(M'^) for all x, the right-hand 
side of (3.5) exists for all z G (C^)c and x G .D. Hence, u{x, z) is defined for 
such z, X. 
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Further, we see for all z e (C^)c that 

< const, Iz"*^! 
(3.6) ' 

1/2 / \ 1/2 

\2N)-2^^ I 



< const. |^^|z^M2(2N)2^fc^ 

< const. -R ■ (^{2ky^^ 



1/2 

< oo 



for all z G K2{R)- Besides this (3.5) shows that u{x^z) is analytical in z. 
Thus, we conclude by the characterization theorem (Theorem 3.8) that there 
exists a function U :D ^ such that ?{U{x, z) = u{x, z). Next we want 

to verify the assumptions of Lemma 3.1 for X = U and F = — It is known 
from the general theory of deterministic elliptic PDE's (see, e.g., [7]) that 
for all open and relatively compact in D there exists a C such that 

(3.7) ||u(-,z)||c2+«(y) <C(||Au(-,z)||c<.(y) + \\u{-,z)\\c(y)) 

for all z G (C^)c. Since An = —7i^ and u are bounded on D x K2{R), it 
follows that -^^u^Xjz) is bounded for such x, z. Thus, by Lemma 3.1 U is 

i 

a solution of system (3.1). 

Further, we follow from Lemma 3.18 in [10] that the Bochner integral 
J^dG{x,y)v{x) dx exists in (5)* (see Definition 3.16 in [10]) and that 

(3.8) / G{x,y)iiiy)dy = mJ2 I G{x,y)Ck{y) dy K.^k.i) . 

Then one realizes that the right-hand side of (3.5) is the Hermite transform 
of (3.8). 

So we obtain the following result. 

Theorem 3.2. There exists a unique stochastic distribution process U :D 
(5)*, solving system (3.1). The solution is twice continuously dijjerentiable 
in (S)* and takes the form 

(3.9) U{x)= f G{x,y)v{y)dy = mJ2 f G{x,y)a{y)dyK.,,,r,, 



where m = ||z||l2(,^). 
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We conclude this section with a remark about an alternative approach to 
SPDEs driven by Levy space-time white noise. 

Remark 3.3. Let us briefly describe how the concepts in [28] can be 
used to establish a framework similar to Section 2. Instead of the spaces 
(5)_p, consider the distribution spaces in [28] and instead of the 7^-transform, 
use the 5-transform in [28]. The 5-transform, is of the form 

S{F){cP) = {{F{u),e{(t>,oo))) 

for distributions F and for (p in an open neighborhood of zero in S{X), 
where the function e{(j),uj) is as in (2.11) and where ((•, •)) is an extension of 
the innerproduct on L^{fi). Moreover, the process ^(x) can be replaced by 

V{x) := {Ci{u}),z5a:) 
and the white noise A'' can be defined by 

iV(x,z):=(CiH,(5(,,,)), 

where 5y is the Dirac measure in a point y. Further, by the properties of 
the 5-transform (see [28]) one can prove a similar result as Lemma 3.1. 
Moreover, the 5-transform of V{x) := {Ci{u)), z6x) is 

5(r? (x) )(</.) = / cPix,zMdz) 

(see proof of Proposition 7.5 in [28]). Hence, we can solve system (3.1) by 
finding a function u such that 

Au{x,(l)) = — / (f>{x, z)iy{dz), x^D, 
JRo 

u{x, (/>) = 0, X G dD. 

The obvious candidate for u is given by the Green function G: 

u{x,4>)= / G{x,y) / (j){x,z)u{dz) dy. 
Jwi JRq 

Hence, the solution is given by the inverse tS-transform, yielding the same 
result as in Theorem 3.2 for all Levy measures. Moreover, within a similar 
setting one can solve more general versions of the problem. However, the 
use of the 7^-transform has some advantages. For instance, it enables the 
application of methods of complex analysis. 
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4. Discussion of the solution. As mentioned in the Introduction, our 
interpretation of the solution U{x) = U{x,uj) of (3.1) is the following: 

For each x we have U{x,-) G (5)* and x i-^ U{x) satisfies (3.1) in the 
strong sense as an (5)*-valued function. 

Regarding the solution itself, given by (3.9), the interpretation is the 
following: 

For each x, U{x) is a stochastic distribution whose action on a stochastic 
test function / e (5) is 

(4.1) {U{x),f)= I G{x,y){v{y),f)dy, 
where 

^ k>i ' 

= mY,^kiy)E[K^.,,.^J] [see (2.40)]. 
fe>i 

Relation (4.1) gives rise to the intepretation that the solution U{x,ij) takes 
averages for all x. 

In general, we are not able to represent this stochastic distribution as a 
classical random variable U{x,uj). However, if the space dimension d is low, 
we can say more: 

Corollary 4.1. Suppose d <3. Then the solution U{x, •) given by (3.9) 
in Theorem 3.2 belongs to L'^{ij) for all x and is continuous in x. 
Moreover, 

(4.2) U{x)= f G{x,y)dr]{y). 

Proof. Since the singularity of G{x,y) at y = x has the order of mag- 
nitude \x — yP""^ for d > 3 and log i ^ for d = 2 (with no singularity for 
d = l), we see by using polar coordinates that 

/ G^{x,y)dy<C f\^^^~'^\'^~^dr 
Jd Jo 

= [\^-'^dr 
Jo 

< oo 

for d<3. Hence, by Remark 2.15 we get 

U{x)= I G{x,y)oViy)dy= [ G{x,y)dr]{y) 
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and 

r2 



E[U\x)]=E 
for d < 3. □ 



G{x,y)dr]{y) 



M f G^{x,y)dy <oo 
Jd 



Some remaining natural questions are the following: 

Ql. If (i > 4, is U{x) G LP(//) for some p = p{d) > 1? 
Q2. Is it possible to prove that equation (3.3) also has a solution U{x) = 
U{x,uj) of Walsh type, that is, such that, for some n, 

xi-^U{x,uj) e H~'^{D) for a.a. 6J 

and U{x) solves (3.1) in the (classical) sense of distributions, for a. a. a;? [See 
(1.3) in the analogous Brownian motion case.] 
We will not pursue any of these questions here. 
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